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ABSTRACT 

The source HESS J1813-178 was detected in the survey of the inner Galaxy in TeV 7-rays, and a composite 
supernova remnant (SNR) G12. 8-0.0 was identified in the radio band to be associated with it. The pulsar 
wind nebula (PWN) embedded in the SNR is powered by an energetic pulsar PSR J1813-1749, which was 
recently discovered. Whether the TeV 7-rays originate from the SNR shell or the PWN is uncertain now. 
We investigate theoretically the multiwavelength nonthermal radiation from the composite SNR G12. 8-0.0. 
The emission from the particles accelerated in the SNR shell is calculated based on a semianalytical method 
to the nonlinear diffusive shock acceleration mechanism. In the model, the magnetic field is self-generated 
via resonant streaming instability, and the dynamical reaction of the field on the shock is taken into account. 
Based on a model which couples the dynamical and radiative evolution of a PWN in a non-radiative SNR, the 
dynamics and the multi-band emission of the PWN are investigated. The particles are injected with a spectrum 
of a relativistic Maxwellian plus a power law high-energy tail with an index of -2.5. Our results indicate that 
the radio emission from the shell can be well reproduced as synchrotron radiation of the electrons accelerated 
by the SNR shock; with an ISM number density of 1 .4 cm~^ for the remnant, the 7-ray emission from the 
SNR shell is insignificant, and the observed X-rays and very high energy (VHE) 7-rays from the source are 
consistent with the emission produced by electrons/positrons injected in the PWN via synchrotron radiation 
and inverse Compton (IC) scattering, respectively; the resulting 7-ray flux for the shell is comparable to the 
detected one only with a relatively larger density of about 2.8 cm"^. The VHE 7-rays of HESS J18 13-178 can 
be naturally explained to mainly originate from the nebula although the contribution of the SNR shell becomes 
significant with a denser ambient medium. 
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1. INTRODUCTION 

The VHE source HESS J1813-178 was discovered with the 
High Energy Stereoscopic Sys tem (HESS) in a survey of the 
inner Galaxy in VHE 7-rays (jAharonian et al.ll2005al) . The 
VHE 7-ray image obtained with the HESS shows a pointlike 
source with an extension of ~ 2.2', and the ob served spectrum 
has a hard photon index 2.09 ± 0.08 (.Aharonian et al.ll2006l) . 
HESS J1813-178 was also detected in VHE 7-rays with the 
Major Atmospheric Gamma Imaging Cerenkov (MAGIC) 
telescope. The differential flux given by MAGIC between 0.4 
and 10 TeV c an be well describe d by a power law with a index 
of -2.1 ±0.2 (I Albert et al.l2006l) . which is consistent with the 
result obtained with the HESS. 

Firstly, HESS J1813-178 was unidentified, and it was as- 
sumed to be a "dark particle accelerator" since no counterpart 
at lower frequencies was reported ever However, the shell 
of the SNR G12.8-0.0 associated with HESS J1813-178 was 
discovered with a diameter of ~ 2.5' in a new low-frequenc y 
Very Large Array (VLA) 90 cm survey (iBrogan et al.ll2005h . 
The non-thermal radio flux densities of the shell-type SNR 
are 0. 65 ±0.10 and 1.2 d= 0.08 Jy at 20 and 90 cm, respec- 
tively dBrogan et al.ll200a) . A highly absorbed X-ray source 
AX Jl 8 13- 178, for which the column density is 10^"^ cm"^, 
detected with ASCA is spatially coincident with the SNR. The 
X-ray emission extending to 10 keV with a sharp cutoff be- 
low 2 keV is primarily nonthermal, which can originate either 
from the SNR shel l or form the pulsa r wind nebula (PWN) 
inside the remnant (iBrogan et al.ll2005h . Due to the high col- 
umn density derived from the ASCA data, a distance of > 4 
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kpc is derived for the source (Brog an et al.ll2005h . 

Moreover, a soft 7-ray source, IGR J18135- 1751, was dis- 
cover ed as the counterpart of HESS J1813-178 (lUbertini et al.l 
l2005h . It is persistent with a 20-10 keV luminosity of 
5 X 10^"^ erg s"' for a distance of 4 kpc. lUbertini et aH (l2005h 
argued that the observed properties of the source in the radio 
and X-ray bands can be explained with the assumption that 
the source is a pulsar wind nebula embedded in G12. 8-0.0. 
Furthermore, high-angular resolution X-ray observation with 
XMM-Newton shows that the X-ray emitting object appears 
as a compact core located in the center of the radio shell- 
type SNR G12.8-0.0 (Funk et al. 2007). They argued that the 
source is a composite SNR since the central object shows mor- 
phological and spectral resemblance to a PWN. 

The observation with Chandra on the SNR G 12.8-0.0 in- 
dicates the X-ray source is a point surrounded by struc- 
tured diffuse e mission that fills the interior of the radio shell 
(iHelfand et al.l l2007h . The compact source has a spectrum 
characterized by a power law with an index of ~ 1.3, typi- 
cal of young and energetic rotation-powered pulsars, and the 
morphology of the diffuse emission strongly resembles that of 
a pulsar wind nebula (Helfand et al. 2007). Recently, an en- 
ergetic pulsar PSR J1813-1749 with a period of -- 44.7 ms, a 
characteristic age of 3.3 - 7.5 kyr, and a distance of 4.7 kpc by 
assuming the association with an adjacent young stellar clus- 
ter, was discovered in a long, continuous XMM-Newton X- 
ray timing observation (Gotthelf & Halpern 2009). The pul- 
sar was found to be associated with the SNR G12. 8-0.0, and 
it powers the PWN ( Gotthelf & Halpern 2009 ). 

High-energy 7-rays can be produced either from SNR 
shells in which particles are accelerated to relativistic through 
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the first Fe rmi pr ocess (e.g., lAharonian et al.l l2005b l 120071: 
Berezhko & Volkl 120 06: Fan g & Zhang 2008; Fan g et alJ 
2009i: iMorhno. Amato & Blasll2009l) . or from PWNe pow- 



ered by the pulsars insid e them (e.g.. IVolpi et al J 120081: 
IZhang. Chen & Fang 2008t iGelfand. Slane & Zhang! l2009l) . 
Although the source HESS J1813-178 is pointlike in VHE 7- 
rays, the possibility of the VHE 7-rays originating from the 
shell of the remnant cannot be ruled out given the size of the 
SNR, the angular resoluti on of the HESS te lescope, and the 
depth of the observations dAlbert et al.l 12006). In this paper, 
we study the multiband nonthermal emission from the shell 
of the SNR G12. 8-0.0 and the nebula inside it. The emission 
from the particles accelerated by the SNR shock is investi- 
gated based on a semianalytical method to the nonlinear diffu- 
sive shock acceleration mechanism with a free escape bound- 
ary proposed by Caprioli et al . (2010b), in which the ampli- 
fied magnetic field due to resonant streaming instability in- 
duced by cosmic rays and the dynamical feedback of this self- 
generated magnetic field on the shock are taken into account. 
On the other hand, the dynamics and the multi-band radiative 
properties of the P WN are investiga ted basically according to 
the model in Gelfan d. Slane & Zhan g (2009), which can self- 
consistently describe the dynamical and radiative evolution of 
a pulsar wind nebula in a non-radiative supernova remnant. 
Recently, based on the long-term two dimensional particle- 
in-cell simulations, Spitkovsky (2008) found that the parti- 
cle spectrum downstream of a relativistic shock consists of 
two components: a relativistic Maxwellian and a power law 
high-energy tail with an index of -2.4 ± 0. 1 . Different from 
iGelfand. Slane & ZhangI (12009 ^. in which a single power law 
injection spectrum for the electrons/positrons is employed to 
discuss the radiative properties during different phase of the 
PWN, in this paper we argue that the high-energy particles are 
injected with the spectrum of a relativistic Maxwellian plus a 
power-law high-energy tail during the evolution, and a kinetic 
equation is used to obtain the energy distribution of the parti- 
cles. 

2. MODEL AND RESULTS 

In this section we describe the physics of our model for the 
particles accelerat ed b y a shock and its multiband nonthermal 
emission (Section lzTI ). the dynamics and multiwavelength ra- 
diation of a PWN in the nonradiative shell (Section |22]l, and 
its implementation to the SNR G12. 8-0.0 (Sectionl23b. 

2.1. Particles accelerated by the SNR shock and its emission 

The pitch-angle averaged steady-state distribution of the 
protons accelerated at a shock in one dimension satisfies 
the diffusive transport equation (Malkov & Druryi i200 h .Blasi, 
l2002l;lAmato. Blasi & Gabicill2008l) . 
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+ Q{x,p) = Q, (1) 



where the coordinate x is directed along the shock normal 
from downstream towards upstream, D is the diffusion co- 
efficient and u is the fluid velocity in the shock frame, which 
equals M2 downstream {x > 0) and changes continuously up- 
stream, from Ml immediately upstream (x = 0~) of the sub- 
shock to M() at far upstream. For the Bohm diffusion, D = 
pc^/(3eB), where B is the local magnetic field strength. The 
maximum momentum of the protons accelerated by a SNR 



increases with time in the free-expansion phase of the SNR 
due to the efficient magnetic field strength and the constant 
shock speed. After the beginning of the Sedov-Taylor phase, 
the maximum momentum drops with time due to the de- 
crease of the shock speed and the efficiency of the magnetic 
field amplification (C aprioli et al. 2009. 2010 a). In this case, 
particles with higher momentum will escape from the SNR, 
and this phenomenon can be mimicked by imposing a free 
escape boundar y at a location up stream of the shock, i.e., 
f(-xo,p) = (Caprioli et al. 2010a'F). 

With the assumption that the particles are injected at im- 
mediate upstream of the subshock, the source function can be 
written as Q{x,p) = Qq{p)6{x). For monoenergetic injection, 
Qoip) is 

Q0(P)= /_ 2 Oip-pinj) , 
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where pi„j is the injection momentum, Wgas.i is the gas 
density at x = 0"*" and 77 is the fraction of particles in- 
jected in the acceleration process. With the injection 
recipe known as thermal leakage, 77 can be described as 
77 = 4(R,nh-l)E^e<'/3y^^ (iBlasi. Gabici & Vannonil llool 
lAmato. Blasi & Gabicill200 8'). where R^ub = u\/u2 is the com- 
pression factor at the subshock and ^ is a parameter of the 
order of 2-A describing the injection momentum of the ther- 
mal particles in the downstream region (pinj = i^/?th,2)- Px\\2 = 
{Inipk^Ti)^/^ is the thermal peak momentum of the particles 
in the downstream fluid with temperature T2, nip is the pro- 
ton mass and k^ is the Boltzmann constant. The downstream 
temperature T2 is ca lculated with equations (10) and (11) in 
ICapriohetal] (120091) . A relatively large ^ = 3.8-4.1 is usually 
employed to investigate the rad i ative properties of SN Rs (e.g., 
iMorlino. Amato & Blasil [20091: iMorlino et al.ll2009l) . and we 
adopt ^ = 3.8 in this paper. 
The normalized pressure in cosmic rays is 



Pcix) ■■ 



All 
3poul 



dpp v{p)f(x,p), 



(3) 



P'mi 



where po is the gas density far upstream of the shock. Mag- 
netic fields can be generated by streaming instability induced 
by the accelerated particles. With the assumption of the tur- 
bulence is generated by the resonant streaming instabihty, the 
normalized p ressure of the amplified m agnetic field can be 
described as (ICaprioh et al.ll2009ll2010"bh 
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4mo f/3/2(x) ' 
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where U{x) = u(x)/uo, va = -So/(47rpo)'''^ is the Alfven veloc- 
ity, Bo is the background magnetic field strength. Then the 
strength of the amplified field is SB(x) = (Sn poUQPi^.(x)y ^^ , 
and the magnetic field dow nstream of the shock is further 
enhanced by B2 = RsubBi (iMorhno. Amato & Bla"sil [2009), 
where Bi = dB(0) is the amplified magnetic field immediately 
upstream of the subshock. The effect of turbulent heating is 
ignored in this paper because the properties of the damping of 
the magnetic field is still uncertain now, and a new parameter 
is needed ev en with the damping p henomenologically taken 
into account (IC aprioli et al.ll2010ah . Neglecting the effect of 
turbulent heating, the normalized pressure of the background 
gas is 

U(x)-^ 
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Fig. 1. — Particle spectra at the shock and the escape fluxes (p^</'csc/«o) 
with 5 = 3.8, To = 10" K, no = 1-4 cm-^, Bq = 5^lG, hq = 1 x 10** cm s"', 
xq = 0.18 pc (solid line); xq = 0.9 pc (dashed line); and xo = 1.8 pc (dotted 
line). 
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Fig. 2.— Pcix) (solid line), U(x) (dashed line), and 5B|lm^lG (dotted 
line) in the upstream region for.vo = 0.9 pc. Others are the same as Fig[T] 

where Mo is the sonic Mach number of the shock, 7 = 5/3 
is the adiabatic index. The total compression factor 7?tot is re- 
lated with the compr ession factor at the subshock/?sub through 
JCaprioli et al.ll2009l) 
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[l+/?,ub(2/7-l)]. 



(6) 



(7) 



Given a value of ^, a temperature far upstream of the shock 
To, a shock velocity uq, a background magnetic field Bq, and 
xo, the particle spectrum at the shock f\ (p) and the escape flux 
0esc(p) = -\ p(x, p)df/dx]x„ can b e obtained with the method 
proposed in lCaprioli et alj (1201 Obi) , in which the diffusion co- 
efficient is calculated in the self-generated magnetic field in- 
duced by resonant streaming instability of the accelerated par- 
ticles. In order to show how the protons spectrum is affected 
by the choice of xq, we plot in FiglUthe particle spectra at the 
sub-shock position and the escaping fluxes for three different 
values of Xq: 0.18 pc (solid line); 0.9 pc (dashed line); and 
= 1.8 pc (dotted line). The other parameters are afl fixed to: 



( = 3.8, To = 10"* K, «o = 1-4 cm-^ Bq = 5^G, mq = 1 x 10** 
cm s"' . The particle spectrum for each xo is cut off at a max- 
imum momentum pm^x, which increases with the value of xq. 
U(x), Pc{x) and SB in the upstream region for xq = 0.9 pc are 
indicated in Fig|2] The amplified magnetic field immediately 
upstream of the subshockBi is 56 pG, which is significantly 
stronger than the background magnetic field. 

The electrons have the same spectrum of the protons up 
to a maximum energy determined by synchrotron losses. 
Note that the spectrum of the accelerated electrons at the 
shock around the cutoff momentum p,mx.e is difficult to ob- 
tain in a fully nonlinear scenario. With the test particle 
approximation, the energy spectrum of electrons acceler- 
ated by SNR shocks for a strong shock cari be described as 
dZirakashvih & Aharonianl2007l;lBlasi|[2010h 

fe(x,p) = Kepf(x,p) [1 +0.523(/7//?max,f)^] exp(-///?^^^ J. 

(8) 
However, the spectrum cut off by a simple exponential is also 
empl oyed to reproduc e the multiband emission from SNRs 
(e.g.. iFang et alj|2009h . and we use this form in this paper, 
i.e., 

fe{x,p) = Kepf{x,p)exp{-E{p)/E^^„,,), (9) 

where E(p) is the kinetic energy of the electrons, and the elec- 
tron/proton ratio Kep is treated as a parameter The choice of 
the cut-off shape of the spectrum has no influence on our re- 
sults because in this paper we interpret the observed X-ray 
from emission G12. 8-0.0 as due to the PWN inside the rem- 
nant. The maximum energy of electrons results equating the 
synchrotron loss time with the acceleration time. In the con- 
text of non linear shock acceleratio n an approximate solu tion 
is given by the following equation dBerezhko et al.ll2()02l) 
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MeV, 



(10) 

where the magnetic field immediately upstream of the shock, 
B\, is assumed to be compressed by a factor /?siib- 

Assuming the accelerated particles distribute homoge- 
neously and most of the emission is from downstream of the 
shock, and using the distribution function at the shock to rep- 
resent the particle distribution in the whole emitting zone, the 
volume-averaged emissivity for photons produced via p-p in- 
teractions can be written as 



e(£) = 47rnga,s / dEpJp{Ep) 



da(E , E„) 



dE 



(11) 



where Ep is the proton kinetic energy, Wgas = Rtotno is the 
gas number density downstream of the shock, and Jp{Ep) = 
vp^fo(p)dp/dEp is the volume-averaged proton density and v 
is the particles' velocity. We use the different ial cross- s ection 
for photons da(E,Ep)/dE presented in Ka mae et alj ( l2006h 
to calculate the hadronic 7-rays produced via p-p interaction. 
Finally, the photon flux observed at the earth can be obtained 
with 

VQ(E) 
F(E)=£^, (12) 

where d is the distance from the earth to the source and 
V is the average emitting volume of the source. For the 
accelerated protons, the e mitting volume can be estimated 
by Vp « (47r/3)7;;^„//?,o, (lElUson et alJ l2000l) . here /?,„,- is 
the radius of the SNR. For the electrons, the thickness 
of the emitting region can be estimated by solving the 
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diffusion-convection equation downstream of the shock, i.e., 
U2{dfe/dx) = Dd^ Idx^- fe/ Tsvn, where Tsvn is the tinie scal e 
of the synchrotron radiation dMorlino. Amato & Blasill2009l) . 
The solution sh ows fe oc exp(-x/j^rini), and the spatial scale 
Rrim is given by (iBerezhko & Volkll2004l) 



Riimip) = 



2D2/U2 



^\+AD2/u2T, 



syn 



1 



(13) 



As a result of the synchrotron losses, electrons with relatively 
high energy are confined into a thin rim behind the shock and 
the total emitting volume is smaller than Vp. A break in the 
spectrum occurs for p = pi, defined as the momentum where 
the time scale of the losses equals to the age of the remnant, 
i.e. Tsyn(pt) = fsni- In the steady state, the volume of the elec- 
trons with higher energies is smaller than that of the protons 
due to the strong synchrotron losses. Hence Veip) = Vp for 
p < pi, while Ve(p) = VpRnUp)/iRsn,/3Rtot) for p > /?„ where 
Pt is determined by A^TR^^^^.Rn,J^(p) = Vp. 

2.2. Dynamics and radiative properties of a PWN inside a 
nonradiative SNR 

A PWN is powered by the pulsar which dissipates its rota- 
tional energy into the nebula. The spin-down luminosity of a 
pulsar with a rotation period of P ev olves with time as (e.g., 
iGaensler & Slanell2006l: [slanai2008h 



£(f) = £o 1 + 



To 



(14) 



where tq is the spin-down time scale of the star, £0 is the 
initial spin-down power, n is the braking index of the pulsar, 
which is equal to 3 for magnetic dipole spin-down. 

High-energy particles and magnetic fluxes are injected 
into the PWN from the terminate shock located where 
the ram pressure of the unshocked wind is equal to that 
of the nebula. In this paper, we assume the spin- 
down power is distributed between electrons and positrons 
(Es(t) = rieEjt)), and magnetic fie lds (Eb = mEit)) (e.g., 
iGelfand. Slane& Zh ang 2009). In iGelfand. Slane & ZhangI 
(|2009), they used a simple power-law injection spectrum 
for the electrons/positrons to discuss the radiative proper- 
ties during different phase of the PWN. However, a bro- 
ken power-law spectrum is usually needed to reproduce the 
non-thermal emissio n of a PWN with multi-band obser- 
vations (e.g., lAt oyan & Aharonian 1996; Venter & d e Jagerl 
I2006t ISlaneetal.2008: Zhang. Chen & Fa ng 2008). Re- 
cently, based on the long-term two- dimensional particle-in- 
cell simulations, ISpitkovskvl ( 12008!!) found that the particle 
spectrum downstream of a relativistic shock is a Maxwellian 
plus a power-law tail with a index of 2.4 ± 0. 1 , and this spec- 
trum also wa s used to investigate the multiband emission from 
PWNe (e.g.,|Fang^Zhang 2010; Slane et al. 2010). We as- 
sume high-energy particles injected in a PWN are accelerated 
by the termination shock (TS) which is typically relativistic 
with a lorentz factor ^10^ with respect to the pulsar wind up- 
stream of the shock. Therefore, we assume the spectrum of 
the high-energy particles injected in the PWN has the form. 



Q(E,t)-- 



C(f)f exp 
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F < F < F 

(15) 

where, a = 2.4 ±0.1, E^, - 2.6 x 10-'*7ts,6 MeV, 7ts,6 is the 
Lorentz factor of the upstream pulsar wind of the TS in 




Fig. 3. — Upper panel: radius of the SNR (Rsm, dashed line), reverse 
shock (i?is, dotted hne), PWN (Rpwn, solid line), and the location of the pulsar 
(Rpulsai, dash-dotted line). Lower panel: magnetic field strength in the PWN 

(^pwn)- 

units of 10^, Emin = fminEb with /min ^ 7, / is normalized by 
^min/^bexp {-E^in/Eb) = /(^min^b)""- C(t) Can be obtained 
with 



C(t)-- 
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^El+frl 
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2-a 



(16) 



Assuming the particles are homogeneously distributed in the 
PWN, and the energy distribution of these particles in the neb- 
ula evolves as 

dN{E,t) d [^^(£^,)]^g(£^^)^ (17) 



dt 



dE 



where E is the energy-loss rate of the particles with energy 
E. Energy-loss mechanisms include synchrotron radiation, IC 
scattering and adiabatic loss. 

The dynamics of the PWN inside the supernova shell 
is calculated basically following the model presented in 
IGelfand. Slane & ZhangI (l2009l) . The model assumes that 
the progenitor supernova ejects material with mass Mgj and 
energy Esn into the ambient matter with a constant density 
PiSM- Assuming the PWN has no influence on the dynam- 
ics of the forward shock and the reverse shock, the velocity 
and the radius of the forward shock and the reverse shock 
of the surrounding SNR a re calculated with the equations in 
iTruelove & McKea (Il999l) . The pulsar dissipates energy into 
the PWN, which sweeps up the supernova ejecta into a thin 
shell surrounding the nebula, and new particles are injected 
into the PWN at each time step. 

The dynamical evolution of the PWN with the parameters 
for the SNR G12.8-0.0 (see Tabl® is shown in FigO The 



TABLE 1 

Input Parameters for the SNR G12. 8-0.0 



Parameter 


Value 


Parameter 


Value 


d 


4.7 kpc 


Eo 


4.2 X 103** erg s"' 


Esn 


0.4 X 10™ erg 


Tsd 


500 yr 


Mej 


3Mq 


n 


3.0 


"0 


1.4cm-3 


r/B 


1x10-4 


Age 


1200 yr 


£.max 


1000 TeV 


To 


10''K 


Eb 


3 X 10^ MeV 


Ba 


5^0 


I'psr 


120kms-' 


Kep 


1.8 X 10-3 







radiuses of the SNR (7?sm), the reverse shock (Rrs), the PWN 
(Rpvin) and the pulsar are indicated by the dashed, dotted, soHd 
and dash-dotted Hnes, respectively. A velocity of 120 km s"' 
is used in the calculation, which has no influence on the re- 
sulting dynamical structure and radiative property of G12.8- 
0.0 because the young pulsar is safely in the nebula now (see 
Section l23T l. Initially, the pressure of the PWN is much big- 
ger than the pressure of the surrounding supernova ejecta, so 
the PWN expands adiabatically into the cold supernova ejecta. 
The ejecta surrounding the PWN is swept up to a thin shell, 
which is decelerated by ram pressu re since its velocity is big- 
ger than the local sound speed (Gelf and et al.l 2007). The mass 
of the PWN Mpwn increases continuously since the shell ex- 
pands faster than the ambient ejecta. This expansion phase 
ends when the PWN collides with the reverse shock of the 
SNR. After the collision, the pressure inside the nebula Pp„,i 
is much smaller than the pressure of the material around it 
Anr(^pwn)- The vclocity of the PWN decreases greatly, and 
finally the PWN is compressed. During this process of com- 
pression, the magnetic field strength in the nebula increases 
significantly. Furthermore, the radius of the PWN decrease 
significantly, and the PWN will expand again when the inner 
pressure eventually becomes bigger than that of the surround- 
ing ejecta. The nebula experiences a series of contractions and 
re-expansions until the SNR enters the radiative phase of its 
evolution. The pulsar moving in the space will leave the PWN 
during the compression, and it can re-enter the nebula when 
the nebula expands again. With the parameters in Tabled] the 
pulsar firstly leaves the nebula at ~ 13700 yr, and re-enters it 
at- 19000 yr. 

2.3. Application to Gl 2.8-0.0 and discussion 

The shell with a diameter of r^ 2.5', corresponding to a ra- 
dius of 1.7 pc for a distance of 4.7 kpc, of the SNR G12.8-0.0 
had been revealed wit h the radio o bservations with a spec- 
tral index of — -0.48 ("Broga n et al.ii2005 ). Moreover, a pul- 
sar wind nebula with a diameter of ~ 80" embedded in the 
SNR was been disclosed in the observations from the X-ray 
to soft 7-ray bands dUbertini et alj 120051: iFunk et all 120071: 
iHelfand et al.1 120071) . An energetic pulsar with a spin-down 
power bigger than 10^^ erg s"' is argued to exist in the nebula 
based on the ASCA and Chandra observations (Helfand et al] 
l2007h . Recently, the energetic pulsar PSR J1813-1749 associ- 
ated the PWN was discovered in the long, continuous XMM- 
Newton X-ray timing observation, and the pulsar has a period 
of ~ 44.7 ms, a spin-down age of 3.3 - 7.5 kyr, and a spi n 
power of 6.8 ± 2.7 X lO" erg s"' (iGotthelf & Halpernll2009h . 

The small radius places the SNR G12. 8-0.0 amongst the 
smallest known SNRs, which also suggests a young age for 
the remnant. In the model, the radius of the SNR is deter- 
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Fig. 4. — Upper panel: the spectral energy distribution for the PWN in 
G12. 8-0.0 at an age of 1200 yr. Synchrotron radiation (dashed line), and IC 
scattering on the CMB (dotted line), IR (dash-dotted line), star light (dash- 
dot-dotted line) and the synchrotron photons (short dashed line) are shown 
in the figure. The solid line represents the whole IC scattering on all the 
soft photons. Lower panel: the spectral energy distribution of the emission 
from the shell of the SNR G12. 8-0.0. Synchrotron radiation (dashed line), 
bremsstrahlung (dotted line) IC scattering (dash-dotted Une) of the acceler- 
ated electrons and p-p interactions (solid line) of the accelerated protons are 
demonstrated in the figure. The radio data are from the VLA, Bonn, Parkes, 
and Nobeyama observatories iBrosan et al. 2005). The X-ray data are from 
XMM-Newton ( Funk et al...2007.) and INTEGRAL ( Ubertini et al._2005.). T he 
HESS flux points in the VHE 7-ray band are from lAharonian et all )2003) . 
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Fig. 5. — A zoomed view of the resulting spectrum and the HESS results 
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mined by the explosion energy, the mass of the ejecta, the 
ambient density, and the age of the system. We find that, with 
Esn = 0.4 X 10^° erg. Me, = 3Mq, and an ambient ISM den- 
sity «() =1-4 cm"-', the resulting radius of the SNR shell is 
^ 1.8 pc at an age of ^ 1200 yr, which is consistent with 
the radius of the radio shell. With a larger explosion energy 
(Esn > 10^" erg) or a smaller ambient density (no < 0.1 cm" ), 
the age of the system, which is constrained by the radius of 
the SNR, is much smaller the radius of the radio shell, and 
the resulting flux of X-ray emission from the PWN is usu- 
ally much higher than the observ ed one by XM M-Newton 
(iFunk et al.ll2007l) and INTEGRAL (iBrogan et al.ir2005) since 
the magnetic field is stronger for a smaller age (see the lower 
panel of FigOJ. A sub-energetic explosion with Zisn < 10^" 
erg and a low ejecta mass Mej = 3Mq was also argued to ex- 
ist from the mult iwavelength study of the newly discovered 
SNR G310.0-1.6 (iRenaud et al.L2009.) . With m = 500 yr and 
E() = 4.2 X 10^^ erg s"', the spin-down power is 3.6 x 10^^ erg 
s"' at an age of 1200 yr from EqlT4l which is a little smaller 
than the present observed v alue, i.e., 6.8 ±2.7 x 10^^ erg s"' 
dGotthelf & Halpernll2009b . The radius of the PWN at this 
age is 0.7 pc, and the ratio ^p„n/^snr is ~ 0.4. The result- 
ing multiband radiation from the PWN in the composite SNR 
G12. 8-0.0 at an age of 1200 yr is shown in the upper panel 
of Fig m with the parameters in Table[T] In the figure we also 
plot the observed data from radio to Te V band for comparison. 
The data sources are reported in the figure caption. For each 
panel, the interstellar soft photons scattered in the IC process 
contain the CMB, IR (Tir = 35 K, Uir = 1 .0 eV cm"^) and star 
light (Ts, = 3000 K, f/st = 1 .5 eV cm"^). The observed emission 
detected with XMM-Newton and INTEGRAL from the X-ray 
to hard X-ray bands can be explained as the synchrotron ra- 
diation from the electrons injected in the PWN, although the 
resulting flux is a little higher than that observed with XMM- 
Newton. In the VHE 7-ray band, the photons are mainly pro- 
duced via IC scattering of the electrons in the nebula on the 
CMB and IR photons from interstellar dust (Fig|5l). 

The upstream temperature for a SNR varies between ^ 
10^ K for a typical ISM up to ~ 10^ K if the SNR 
expands i n the hot bubble generated by the progeni- 
tor's w ind (Chevalier & Liang 1989^. iMorlino. Am ato & Blasil 
l2009t [Berezhko & Volk.2010.) . We assume the SNR G12.8- 
0.0 is expanding in the typical ISM with a temperature of 



^ 10"* K to investigate the multiband nonthermal emission 
from the SNR shell. At an age of 1200 yr, the velocity of the 
SNR shock is ~ 1000 km s"', corresponding to a Mach num- 
ber of ~ 86 for T() = lO'* K. With these parameters, /?tot and 
Rsuh are about 8.1 and 3.55, respectively, and the downstream 
temperature T2 is 3.8 x 10^ K. Moreover, the maximum ener- 
gies of the accelerated protons and electrons are about 60 TeV 
and 4 TeV, respectively. The energy contained in the protons 
accelerated by the shock is 5 x lO'*'^ erg, and it in the electrons 
is 6 X 10"*^ erg. Therefore, about 12% of the explosion energy 
/isn has been converted to the kinetic energy of the particles. 
The amplified magnetic field strength immediately upstream 
of the subshock is 62 /^G, and then the downstream magnetic 
field strength is 220 /iG. With this downstream magnetic field 
strength and K^p = 1 .8 x 10"^, the flux points of the radio shell 
can be reproduced as synchrotron radiation of the electrons 
accelerated by the SNR forward shock wave; whereas the flux 
from p-p collisions of the accelerated protons on the ambient 
matter is significantly smaller than that observed with HESS. 

The ISM number density «(> is an important parameter 
which can greatly influence the multiband nonthermal emis- 
sion from the SNR shell. With a larger no, both the density 
of the accelerated particles and that downstream of the shock 
increase accordingly. As a result, the flux of the p-p colli- 
sions is enhanced with oc «q, and then a smaller K^p is usually 
needed to reproduce the observed radio fluxes for a denser 
medium. In order to reproduce the HESS flux in the TeV band 
via p-p collisions from protons accelerated by the SNR shell, 
a denser medium with a density of 2.8 cm"-* must be used in 
the model (see Fig|6]l. For no = 2.8 cm"^ and the other param- 
eters the same as FigH] the multiband nonthermal emission 
from the PWN is nearly the same as that for «o =1.4 cm"-', 
whereas the 7-ray flux for the p-p collisions in the SNR shell 
is comparable to the HESS result; the downstream magnetic 
field strength is 27 1 /iG, and then a smaller K^p = 4.5 x lO"'* is 
needed to reproduced the observed fluxes in the radio band. In 
this SNR-dominated scenario, the 7-ray photon index of the 
resulting emission is < 2 up to several tens of GeV, whereas 
it is > 2 at higher energies (Fig|6l). On the other hand, the 
index in the PWN-dominated case is < 2 from 0. 1 GeV to 1 
TeV (Fig|4]l. The resulting TeV 7-ray spectral energy distri- 
butions in the two scenarios are all consistent with the HESS 
flux points, so we cannot determine which case is preferred 
by comparison with the detected fluxes in the TeV range now. 
However, the flux in the GeV band in the PWN-dominated 
case is several times higher than in the SNR-dominated case, 
so future detections in this band can give constraints on the 
origin ofthe7-rays from HESS J1813-178. 

Even if a denser medium with «o ^ 2.8 cm"^ is used in the 
model to enhance the p-p collisions in the SNR shell, the IC 
scattering of the electrons/positrons in the PWN is still promi- 
nent enough to contribute significantly to the high-energy 7- 
ray emission from the composite SNR. In the model, the PWN 
evolves in the remnant with an age of ^ 1200 yr, and high- 
energy electrons/positrons are injected continuously into the 
nebula from the energetic pulsar during the evolution. In fact, 
if the emission detected with XMM-Newton and INTEGRAL 
is the synchrotron radiation from the electrons/positrons in- 
jected in the PWN, the detected VHE 7-rays can be easily 
explained by the IC scattering of these electrons/positrons off 
the ambient soft photons. Therefore, it is natural to argue that 
the observed high-energy photons from the X-ray to VHE 7- 
ray bands originate mainly from the PWN inside the compos- 
ite remnant although the contribution of the SNR shell to the 



7-rays becomes significant for a denser medium with a den- 
sity of 2.8 cm""*. 

3. SUMMARY AND CONCLUSION 

The origin of the VHE 7-rays from HESS J1813-178 is in- 
vestigated in this paper Ahhough the source is pointHke in 
VHE 7-rays, a shell origin of the VHE 7-rays cannot be ruled 
out because the PSF of the HES S telesc ope is very close to 
the size of the SNR's shell (e.g., Albert et al j 2006). A SNR 
with a diameter of 2.5' was identified in the radio bands to 
be associated with the VHE 7-ray source. Moreover, X-ray 
observations showed a P WN powered by an energetic pulsar 
is embedded in the SNR (Brogan et al. 2005] lUbertini et"an 
|2005; Funketal. 2007; H elfand et al. 2007) , and the pulsar 
was recently discovered i n a long, continuous XMM -A^ewfow 
X-ray timing observation dGotthelf & Halpernll200 9). We ap- 
ply a model which can self-consistently calculate the multi- 
wavelength nonthermal emission both from the SNR shell and 
from the PWN embedded in the remnant. In the model, elec- 
trons/positrons with a spectrum of a relativistic Maxwellian 
plus a power-law high-energy tail are injected in the nebula 
during the PWN during its evolution inside the SNR; pro- 
tons and electrons are accelerated by the SNR shock wave, 
and the spectrum of the accelerated particles are calculated 
with a semi-analytical non-linear model. Our results indicate 
that: with the parameters in Tabl^D (1) the observed emis- 
sion of the shell in the radio bands can be well explained as 
synchrotron radiation of the electrons accelerated by the SNR 
shock wave, whereas the flux of p-p collisions of the accel- 
erated protons is significantly smaller than the TeV flux ob- 
served with HESS; (2) the observed emission in the X-ray to 
hard X-ray band detected with XMM-Newton and INTEGRAL, 
respectively, can be explained as the synchrotron radiation of 
the electrons/positrons injected in the PWN; (3) the VHE 7- 
rays in the TeV band are mainly produced via IC scattering of 
the electrons/positrons injected in the nebula on the CMB and 



interstellar IR photons. 

With a sub-energetic explosion {E^^ ~ 0.4 x 10^° erg), an 
age of ^ 1200 yr, and an ambient density of 1.4 cm"-', the ra- 
dius of the composite SNR G12. 8-0.0, 1.7 pc, as well as the 
multiband observed nonthermal fluxes for the remnant can be 
reproduced within this model described in this paper The 
present value of the SNR shock's velocity is about 1000 km 
s~', and the VHE 7-rays are predominately produced from 
the nebula inside the remnant via IC scattering. Of course, 
p-p interaction from the SNR shell can be enhanced with a 
denser medium around the remnant (Fig|6ll. Even if this con- 
dition were satisfied, the 7-rays produced via IC scattering 
in the PWN remain significant, hence, in this case, both the 
SNR shell and the PWN would contribute to the production 
of the observed 7-ray flux. In GeV 7-rays, the spectral prop- 
erty of the resulting emission in this SNR-dominated scenario 
(Fig|6ll differs significantly from that in the PWN-dominated 
case (Figgli. The VHE 7-rays of HESS J1813-178 detected 
with HESS can be naturally explained as the IC scattering of 
the electrons/positrons injected into the PWN although the p- 
p collisions become important with a denser ambient medium 
for the remnant. Our study give more insights on the nature 
of the multiband nonthermal emission of the composite SNR 
G12. 8-0.0, even though some assumptions are made in the 
model. 
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